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Using the Kaluza-Klein structure of stationary spacetimes, a framework for analyzing stationary 
perturbations of static Einstein-Yang-Mills configurations with bosonic matter helds is presented. It 
is shown that the perturbations giving rise to non-vanishing ADM angular momentum are governed 
by a self-adjoint system of equations for a set of gauge invariant scalar amplitudes. The method 
is illustrated for SU(2) gauge helds, coupled to a Higgs doublet or a Higgs triplet. It is argued 
that slowly rotating black holes arise generically in self-gravitating non-Abelian gauge theories with 
bosonic matter, whereas, in general, soliton solutions do not have rotating counterparts. 


I. INTRODUCTION 

In the presence of a (stationary) Killing symmetry, the Einstein-Maxwell (EM) equations reduce to a cr-model 
coupled to three-dimensional gravity ^]. This property is, in fact, shared by a large class of theories with scalar 
fields and Abelian vector fields (see [|| for a classification and |^, [Q for some recent applications and references). If 
spacetime admits an additional (axial) Killing symmetry, then the cr-model structure gives rise to total integrability 
of the field equations, provided that the target space is a symmetric space. This has been known for quite some time 
for the EM system jhj and was recently demonstrated by Gal’tsov for EM-dilaton-axion models . 

Since scalar magnetic potentials fail to exist in non-Abelian gauge theories, the cr-model structure - and, in particu¬ 
lar, the property of integrability - are spoiled for self-gravitating Yang-Mills fields. Moreover, the circularity theorem 
0 (which guarantees that spacetime admits a foliation by two-surfaces orthogonal to the integral trajectories of the 
two Killing fields) does not extend to the Einstein-Yang-Mills (EYM) system 1^ (see also |^). The familiar Papa- 
petrou metric does, therefore, not take account of all stationary and axisymmetric degrees of freedom of the EYM 
equations. 

In view of these problems, an analytic approach to the full EYM equations with two Killing fields is likely to be 
extremely difficult. Motivated by recent work of Straumann and Volkov ||ll|| , we pursue a more modest aim in this 
paper, that is, we consider stationary deviations of static EYM configurations with bosonic matter fields. For the 
pure SU(2) EYM system, Straumann and Volkov were able to reduce the relevant perturbation equations to a 
three-dimensional set. In this paper, we present a systematic investigation of stationary perturbations, which reveals 
that the decoupling of a specific set of perturbation amplitudes is a general feature of a large class of bosonic matter 
fields coupled to the EYM system with an arbitrary gauge group. We argue that stationary perturbations are most 
appropriately handled by means of a “3 -I- 1” - rather than a “2 -|- 2” - decomposition of spacetime. (This does, 
in particular, avoid the circularity issue, since the metric is not required to be axially symmetric in the first place.) 
Hence, we use the Kaluza-Klein (KK) structure of a stationary spacetime to analyze arbitrary stationary perturbations 
of static configurations. Within this approach, the non-static deviations are encoded in the KK connection, which is 
related to (the dual of) the twist of the stationary Killing field. 

The KK reduction of the Einstein-Hilbert action yields a three-dimensional gravitational theory coupled to the 
KK scalar field and the the KK connection [^. The latter is described by a gauge potential, which enters the 
effective action only quadratically and only via the field strength. Using a suitable KK decomposition of the YM 
gauge potential, these features are preserved if gauge fields and additional bosonic matter are coupled to gravity. 
More precisely, it turns out that the linear terms in the KK connection enter the reduced action via a (non-minimal) 
conpling to the electric components of the YM field. These observations imply the following two conclusions: (i) 
The stationary EYM equations (coupled to bosonic fields) admit a generalized scalar twist potential, and (ii) the 
non-static, non-magnetic deviations of a static, purely magnetic solution to the EYM equations form a consistent 
subset of all stationary perturbations. Moreover, it is exactly this subset of perturbations, henceforth called purely 
stationary perturbations, which gives rise to a non-vanishing ADM angular momentum. 
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By virtue of the crucial features (i) and (ii), the relevant perturbations (as far as angular momentum is concerned) 
of a static, purely magnetic EYM-Higgs configuration form a formally self-adjoint system for a set of gauge invariant 
scalar amplitudes. For a spherically symmetric SU(2) background, these amplitudes, consisting of the generalized 
twist potential and the (Lie algebra valued) electric YM potential, can be expanded in terms of “isospin” harmonics, 
Cj^. Since only j = 1 contributes to the ADM angular momentum, one finally obtains a standard Sturm-Liouville 
problem for three radial functions. For the twist channel one has j = i = 1, whereas the orbital angular momenta in 
the two YM channels are £ = 0 and £ = 2. The Higgs fields enter the perturbation equations only via a background 
potential, which gives mass to either one (triplet) or both (doublet) YM perturbations. 

For a stationary background, the horizon is a regular singular point of the perturbation equations, which admit 
four acceptable solutions, whereas the corresponding number is three in the asymptotic regime. The fact that the 
perturbation equations admit a six-dimensional fundamental system then yields the conclusion that slowly rotating 
black hole solutions to the EYM-Higgs equations do exist. The corresponding solutions for the pure EYM system were 
recently discovered by Volkov and Straumann [|ll[, who also argued that these configurations cannot be electrically 
neutral. The perturbation equations show that the coupling of isospin and orbital momentum, which is responsible for 
the “charging up” due to rotation, does not need to be effective if bosonic matter is coupled to the EYM equations. 

For solitonic background solutions the origin is a regular singular point of the perturbation equations. The number 
of physically acceptable modes at the center is, however, not sufficiently large to allow for “generic” rotational degrees 
of freedom of self-gravitating bosonic matter coupled to non-Abelian gauge fields. 


II. KALUZA-KLEIN REDUCTION 

We consider the action for self-gravitating non-Abelian gauge fields coupled to bosonic matter, 

S = “YgttG J 

where k = SirG/g'^, G is Newton’s constant, and g is the gauge coupling. The four-forms Cq and £ym are the 
Einstein-Hilbert and the YM Lagrangians, respectively, 

Ca = , Cym = 2Tr{F(^) A . (2) 

Here, and denote the Ricci scalar and the Hodge dual with respect to the spacetime metric The 

one-form is the Lie algebra valued YM gauge potential with field strength A A^^'>. For the 

bosonic matter we shall, for instance, consider a Higgs field H [with potential P{H)] which transforms according to 
some representation U of the gauge group, = dH Ui,{A^^'>)H. In particular, 

Cb = -2Tr|(L>(4)iLy A (3) 

for a Higgs doublet or a triplet in matrix representation [see also Eq. (^ )]. 

Our first aim is to perform the KK reduction of the above action ^ ). At least locally, a stationary spacetime 
[with Killing field dt and corresponding one-form k = —a{dt a)] has the structure M x Y, and admits a 
metric of KK type, 

= —a (dt -I- a) 0 {dt -I- a) -f a~^g . (4) 

Here, a and a are, respectively, a scalar field and a one-form on the three-dimensional Riemannian space (E,g). 
Under coordinate transformations the one-form a transforms like an Abelian gauge potential. The corresponding field 
strength, da, is proportional to the dual of the twist one-form, w = ^ (k A dk) = —^cr^ * da. (Here and in the 
following * denotes the Hodge dual with respect to the Riemannian metric g.) The canonical decomposition of the 
gauge field in terms of a stationary function (j) and a stationary one-form A (both Lie algebra-valued) on (S, g) 
is 


= (j){dt -\- a) A. (5) 

In the following it will be crucial that A^'^'> is decomposed with respect to the orthonormal tetrad field 0° = ^/a{dt a) 
(rather than ^/adt). The reduction of the Einstein-Hilbert action with respect to the stationary metric (^ gives 
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f Cq = J{dt A Lg), where the three-form Lq is the Lagrangian for the KK scalar field a and the Abelian gauge field 
a, effectively coupled to 3-dimensional gravity. Up to an exact differential, one finds 


Lg 


^ da A *da da A *da . 
2a‘‘ 2 


( 6 ) 


The dimensional reduction of the YM action yields an effective YM-Higgs theory, with effective Higgs field 4> and YM 
field strength F = dA + A A A. With f Cym = J {dt A Lym) one has 


Lym = 2Tr 


-|cr(F + (j) da) A -^{F (j) da) 
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Z3(/) A *D(f) 


( 7 ) 


where D denotes the (gauge) covariant exterior derivative with respect to the one-form A on S. Introducing a field 
strength vector with components da and U, the above formulas imply that the stationary EYM system reduces to a 
three-dimensional EYM theory which is non-minimally coupled to a two-component vector of scalar fields (comprising 
combinations of the KK scalar a and the YM scalar (j)). Finally, the evaluation of the Higgs action with respect to 
the gauge potential (||) results in an additional potential term, involving the coupling between the actual Higgs field 
H and the effective Higgs field (j): 


Lb = -2Ti A *DH - ^ (C/*(<())iJ)^ A *t7*(<())7j| - i * P[H]. (8) 

The vacuum Einstein equations are obtained from variations of f Lq with respect to g, a and a. Since Lq is a 
quadratic expression in terms of da, both the effective three-dimensional Einstein equation for g and the equation for 
a contain no linear terms in da. In the presence of YM and Higgs fields this property generalizes in the sense that the 
effective action continues to be quadratic in combinations of da and (j). Hence, the only equations which contain linear 
terms in da and/or (j) are those which are obtained from variations of the effective action, [Lg -\- k {Lym + Lb)], 
with respect to these quantities: 


d * da -I- 4k ct Tr {</> {F + (f>da)}] = 0 , (9) 

D*[a~^D(j)] + a da A *{F + (j)da) = a~'^ * Jb{4>) , (10) 

where Jb{4’) is the bosonic current (zero-form). In particular, one has 

Jb(</)) = - [id, [id,/)]] , and Jb(/)) = ^ (./idtid + idtid/)) (11) 

for a Higgs triplet and a Higgs doublet (in matrix representation), respectively, provided that the latter transforms 
by left multiplication under the action of SU(2). 

Equation m is the electric part of the YM equation. The twist Eq. (|^) assumes the form of a differential 
conservation law. This is due to the fact that the connection a is an Abelian gauge field which - for reasons of 
gauge invariance - enters the effective action only via the field strength da. All stationary self-gravitating matter 
models give, therefore, rise to a generalized twist potential, x, say. It is well known that the twist potential for the 
Einstein-Maxwell system involves the electric and the magnetic potential. The above reasoning implies that the twist 
potential continues to exist in the EYM system, although scalar magnetic potentials cease to do so in non-Abelian 
gauge theories. In fact, Eq. (||) implies the existence of a function x, such that 

(l -I- 4Kcr“^Tr {/)^}) da = a~^ * dx — 4Kcr“^Tr {(pF} . (12) 

(It may be worthwhile mentioning that an explicit expression for the twist potential does not exist for a rotating 
boson star. This is a consequence of the fact that the effective action does contain terms in a itself, since the model 
is not stationary in the strict sense and is, therefore, only gauge invariant under a combined transformation involving 
a and the time coordinate.) 
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III. STATIONARY PERTURBATIONS OF STATIC SPACETIMES 


Let us now consider stationary perturbations of a static (i.e., a = 0) EYM configuration. The above reasoning 
implies that the perturbations da and Scj) do not couple to the remaining metric and matter perturbations, provided 
that the static configuration is purely magnetic. (In this case both a and are first order quantities.) The stationary 
perturbations of a static, purely magnetic spacetime therefore fall into two complementary sets, henceforth called static 
perturbations and purely stationary perturbations. The static set involves only perturbations of fields (metric and 
matter) which are already present in the equilibrium configuration. It is obvious that the restriction to perturbations 
of this kind gives rise to a consistent set of first order equations. The purely stationary perturbations involve those 
fields which vanish for static, purely magnetic configurations. It is an interesting consequence of the above KK 
reduction that the purely stationary perturbations form a consistent subset as well, that is, the twist channel and the 
electric channel do not cause perturbations of the remaining fields. 

It is very intuitive (and will be shown below) that it is precisely the set of purely stationary perturbations which 
gives rise to angular momentum. Hence, we shall now focus on these perturbations, that is, we consider 


5g = 0, 6a = 0, 


5A = 0, 


SH = 0 


(13) 


and 


= Sa 


= 5(j). 


(14) 


The arguments presented above imply that the static equations for g, ct, A and H remain unchanged in first order 
perturbation theory. The perturbation equations for 6a and dtp are obtained from Eqs. (^ and ([^), respectively. 
However, it turns out to be more convenient to use the linearized twist potential^ 6x, rather than da itself. The 
perturbation equation for dx is derived from Eq. (11) by linearizing the integrability condition d{da) = 0, whereas 
the perturbation equation for dx is obtained from Eqs. (10) and (|l^. One easily finds (to first order in dx and d(j)) 


j;:d(-^*ddx] + dfiTrlFJ^} ) = 0, 


(15) 


D 


* Ddcj )) + —F A ddx = 4 :kTt {F dip} A *F + — * JB{d(j)) ■ 


(16) 


The above equations for the scalar perturbations dx and dp form a formally self adjoint system. This is manifest 
for the second order differential operators and for the diagonal potential terms on the right hand side of Eq. (H)- 
The two off-diagonal parts on the left hand sides are easily seen to be symmetric as well. Moreover, dx and dp are 
gauge invariant perturbation amplitudes: This is obvious for dx-, since it is obtained from the Abelian field strength 
d{da). The invariance of dp follows from the infinitesimal transformation law 5p ^ dp -\- Ui,{p)df and the fact that 
p vanishes for the background solution. [We recall that under an infinitesimal gauge transformation, df, one has 
dA^'^'> dA^^'i + D^'^\df).] Before we proceed, we shall briefly argue that the angular momentum of a stationary 
spacetime involves only the purely stationary set of perturbations, governed by eqs. d^) and dH). 

Apart from stationarity, no symmetry requirements have been imposed so far. We shall now assume that spacetime 
admits a second, axial Killing field, and compute the Komar expression for the angular momentum jj. Asymptotic 
flatness implies that only the terms which are linear in a contribute to the Komar integral, 


J 


1 

IGttG 



1 

16ttG 



[a A *dp — p A *da] . 


(17) 


Here, = g^^dx^ is the axial Killing one-form and "0 its projection on S. Since in the asymptotic regime cr —> 1 
and sin^ddip, the first integrand in Eq. ( 0 ) becomes equal to minus two times the second one. Hence, the 

angular momentum becomes 

- wb 4 = leb L_ 

Let us now consider arbitrary stationary, axisymmetric perturbations of a static and axisymmetric spacetime. In 
this case, a is a first order quantity, and the expression for ff involves neither perturbations of the 3-metric g nor of 
the KK scalar field a. Hence, only the purely stationary modes, governed by Eqs. ( p^ and ( 0 ) , contribute to the 
angular momentum. 
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IV. MULTIPOLE EXPANSION 


We now restrict ourselves to spherically symmetric background configurations and perform a multipole expansion 
of the relevant first order quantities (which, for simplicity, are assumed to be axisymmetric). In the unperturbed 
spacetime {MxT ,, we use standard Schwarzschild coordinates and parameterize the metric = —a dt^+ (J~^ g 
in the familiar form 

cr = iVS'2 , (J~^g = N-^dr"^ + ^ 

where N and S are functions of the coordinate r. In the “canonical gauge”, the static, spherically symmetric, purely 
magnetic background YM potential assumes the form 

A = [\ — w{r)] idTr , (20) 

where # denotes the Hodge dual with respect to the standard metric on S'^, and Tr, t,j, are the spherical generators 
of SU(2) (normalized such that [r^,T^] = Tr). (See also [Q for a discussion of symmetric gauge fields with a higher 
rank gauge group.) For a static, spherically symmetric Higgs field we have 

= h(r)r,, = i/i(^) 1, (21) 

where, as before, and denote a Higgs field in the adjoint (triplet) and the fundamental (doublet) rep¬ 
resentation of SU(2), respectively. (We recall that the general spherically symmetric ansatz for a Higgs doublet is 
= ^h{r)\ — g(r)Tr, and that the magnetic gauge potential A involves the additional term w{r)dTr. However, in 
the static case, the field equations imply that one may consistently set g{r) = w(r) = 0; see, e.g., p^ .) 

Let us now consider the multipole expansion for the perturbations. We first observe that the perturbations of 
the metric potential a which contribute to the ADM angular momentum belong to the sector with (total) angular 
momentum j = 1. In fact, as Sa is an axisymmetric one-form on the spherically symmetric manifold S, this has an 
expansion of the form 


Sa = [aj* dYj + (3j Yj dr + 7 ^- dYj ] 
j 


( 22 ) 


where the coefficients are functions of the radial coordinate r, and Yj is a short hand for the spherical harmonics Yjo. 
Since the integrand in the Komar expression ( p^ ) is proportional to dYi , the orthogonality of the spherical harmonics 
implies that only the term proportional to * dYi in the expansion for da gives a non-trivial contribution. Hence, as 
claimed, the sector describing infinitesimal rotations consists of the purely stationary perturbations with total angular 
momentum J = 1. 

Next, we evaluate the perturbation equations (lA), (|^) for the background fields A and H (given in eqs. ( ^ ) and 
(il), respectively), which are easily seen to be symmetric under parity. To this end, we first expand the electric YM 
perturbation 5(j) in terms of the “isospin” harmonics which, after suitable identifications, are proportional to the 
standard vector harmonics Y-l,: 

J IIL 


= TA VbY,^ , = T ^ (2j + 1 ± I) A Y,m + 5^^ VeY^m , (23) 

where capital Latin letters refer to indices with respect to the orthonormal frame 9'^ = dd, 9'^ = sind dip on S^. 
The harmonics have parity (—1)^ and are, of course, eigenfunctions of the Laplacian A = *d*d on S'^ with 
eigenvalues —£{£ -1-1). It is not hard to see that the symmetry under a parity transformation implies that the odd 
parity component of S(j) decouples. Moreover, this does not contribute to the ADM angular momentum, since the 
parity of the corresponding variation of a is also odd [see Eq. ®]- Thus, the axial perturbations which are relevant 
to infinitesimal rotations can be parameterized in terms of three scalar functions x(r), y{r) and z(r): 


Sx = V^x{r) Yi , 


S(j) = y(r) TrYi + z{r) Yi . 


(24) 


At this point, it is a straightforward task to derive the perturbation equations for the vector-valued function 
= {x^ y, z)'^ from Eqs. dH) and dl)- The rotational deviations are governed by the following Sturm-Liouville 
equation: 
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( 25 ) 


+ J + Bd - dB'^ + pjv = 0 , 

where d denotes the differential operator 


a/ ^ ^ i ^, ( 26 ) 

and S is defined in Eq. dil)- The first two terms originate from the differential operators D (a and 

d [a~'^ *d6x), which give rise to the matrix-valued background functions 


/-a-i 

0 

0\ 

1 ( 

' -2cr-i 0 

0 \ 


0 

1 

0 , 

0 2(r(;2 -p 1) 

-2v^n; . 

(27) 

V 0 

0 

1/ 

^ \ 

CNl 

1 

o 

^2 -1- 1 / 



(Note that for w —> 1 and a —!■ 1 the eigenvalues of J become —2, 0, 6, which reflects the fact that the twist channel 
has angular momentum j = ^ = 1, whereas the orbital angular momentum of the YM perturbations is 0 and 2.) For 
the differential coupling between the twist potential and the gauge fields we obtain (in units with k/2 = AmGj = 1) 

/ 0 da~^{vP‘ — 1) 0 \ 

Bd - dB'^ = 2 \ -\)d 0 0 . (28) 

\ 0 0 0 / 

Finally, the potential matrix P is given by 

„ / 0 0 V2w' \ 

P = -- I 0 2 (^2 - 1)^-2 0 ] + Ph, (29) 

^ \ w' 0 2 cr Ul'2 / 


where the background Higgs field enters the perturbation equations only via the matrix Ph, which becomes 


p(3) _ 


a 


0 0 0 \ 

0 0 0 

0 0 /i 2 / 


and 


4cr 


0 0 0 \ 
0 /l 2 0 

0 0 /l 2 / 


(30) 


for a Higgs triplet and a Higgs doublet, respectively. 

In order to discuss the pulsation equations one needs the behavior of the background quantities N, S, w and h. These 
are subject to the static, spherically symmetric EYM-Higgs equations, which are most conveniently obtained from the 
effective Lagrangian. For the gravitational part one finds (up to an exact differential) = AS^dt A dr A dfl, 

where 2m(r) = r[l — iV(r)]; see, e.g. Also evaluating the effective Lagrangians (|^) and (||) [with a = 0 and (/i = 0] 
immediately gives the static, spherically symmetric action (using again k/2 = 1 ) 






(^2 — 1)2 
2^2 


-P{h) + Q{w,h) \ Sdr , 


(31) 


where P{h) denotes the Higgs potential, and the interaction potential Q{w, h) is given by 

Q'^^\w,h) = h? , and Q^‘^\w,h) = (1 — , (32) 

for a Higgs triplet and a Higgs doublet, respectively. Variation of S with respect to m and S yields the relevant 
Einstein equations, whereas variation with respect to w and h gives the magnetic YM-Higgs equations. Using the 
background equations enables one now to analyze the perturbation equations in the vicinity of the origin, the horizon 
and in the asymptotic regime. In the following section we present the results of a systematic discussion. 


V. ROTATING BLACK HOLES 

We start by discussing the behavior of perturbations near the horizon, rn, of a given black hole background. If 
the unperturbed solutions are analytic in a neighborhood of the horizon, then rn is a regular singular point of the 
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perturbation equations. Local properties of the solutions can, therefore, be analyzed by means of standard techniques. 
In particular, the number of physically acceptable solutions is easily determined: The perturbation equations for the 
EYM system coupled to a Higgs doublet or a Higgs triplet admit precisely four independent solutions which are 
admissible near the horizon (provided that the unperturbed black hole is not extreme). 

Next we consider the asymptotic regime, r —>■ oo. Near infinity, the background solutions with a Higgs field 
in the adjoint representation approach the embedded Reissner-Nordstrom solution with magnetic charge = 1: 
w « 0, and \h\ « v, where v is the vacuum expectation value of the Higgs field. Similarly, the unperturbed solutions 
with a Higgs field in the fundamental representation approach the embedded Schwarzschild solution: |r(;| « 1, and 
\h\ « V. (The Abelian nature of the matter fields becomes manifest after a suitable gauge transformation.) It is 
straightforward to verify that the leading asymptotic behavior of the perturbations remains unchanged if a given 
background solution is replaced by its “asymptotic Abelian part”. Within this approximation, the perturbation 
equations simplify considerably in the asymptotic regime: For a Higgs triplet, the “massive” perturbation channel 
decouples, and the remaining two equations have a regular singular point at infinity. For a Higgs doublet, the 
asymptotic system can even be decoupled completely. For both types of Higgs fields it is, therefore, readily verified 
that precisely three independent solutions exist which are physically acceptable near infinity. 

Since the background configurations are continuous for rn < r < oo, the above-defined local solutions have 
extensions with a range of definition containing the whole interval rn < r < oo. By construction, these extensions 
span the subspaces of global solutions which are acceptable near the inner and the outer boundary point, respectively. 
Since these solution-subspaces have dimension three and four, respectively, and since the dimension of the total 
solution-space is six, the intersection of the subspaces is (at least) one-dimensional. Thus, physically acceptable 
global solutions of the perturbation equations always exist for the EYM-Higgs system. 


VI. ROTATING SOLITONS 

Like in the black hole case, the perturbation equations for soliton background solutions have a regular singular 
point at the inner boundary point, r = 0, provided that the unperturbed solutions are analytic in a neighborhood of 
the origin. In the vicinity of this point, the leading behavior of perturbations is completely fixed by the “centrifugal 
barrier”, Jjr"^. It is, therefore, straightforward to verify that precisely three independent solutions exist which 
are globally defined and physically acceptable near the origin. In the asymptotic regime, r ^ oo, the behavior of 
perturbations is the same as in the black hole case. Hence, the global solutions of the perturbation equation which are 
admissible near both boundary points are given by the intersection of two solution-subspaces, each of which is three- 
dimensional. Since the intersection of two three-dimensional subspaces of a six-dimensional linear space generically is 
trivial, we are led to the conclusion that soliton solutions of the EYM-Higgs system generically do not admit rotational 
excitations. 


VII. CONCLUDING REMARKS 

Both, the general structure and the main features of the perturbation equations are dominated by the EYM part 
of the system. It is, therefore, natural to expect that the above results, derived for the SU(2) EYM-Higgs system, 
continue to hold for a class of EYM systems with higher rank gauge groups, and more general bosonic matter 
fields. Hence, we conjecture that bosonic EYM black holes always have rotating counterparts, whereas bosonic EYM 
solitons generically do not admit infinitesimal rotations. The approach presented in this paper offers the possibility 
for a systematic study of these conjectures, which, in case they should turn out to be correct, raise the important 
question about the physical mechanism preventing bosonic solitons from rotating. 
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